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An unified approximated solution for symmetric Skyrmions was proposed for the SU(2) Skyrme model for baryon numbers up to 8, 
which take the hybrid form of a kink-like solution and that given by the instanton method. The Skyrmion profiles are examined by 
computing lowest soliton energy as well as the static properties of nucleons within the framework of collective quantization, with a 
good agreement with the exact numeric results. The comparisons with the previous computations as well as the experimental data 
are also given. 

Keywords: Skyrme Model, Soliton, Nucleons 
PACS: 12.38.-t, 11.15.Tk, 12.38.Aw 



Oh 
i 

Oh 
<D 
43 



> 

m 

cn 

On 
O 
O 



X 



1. Introduction 

Due to its substantial connections with the baryon phe- 
nomenology, the Skyrme model|Q]], which appears as the ef- 
fective field theory of QCD at low-energy limit, has been vis- 
ited extensively since 1980's (see, |2|], for a review), of which 
the soliton solution, known as Skyrmion, has been often used 
to explore the baryon and barynic system spectrum and their 
static properties, particularly, in the case of the small baryon 
numbers (B) Owing to the high nonlinearity, 

the most of Skyrmion studies was confined to the numerical 
approach JH [Io|, [l, 0, Q] . It is still worthwhile, however, to 
seek the approximated analytic soluitions of Skyrmions consid- 
ering its extensive applications in baryon phenomenology and 
nucler physics. The some of efforts have been made in find- 
ing the approximated solution to the massless Skyrme model 
lflu l^ [l2lll3ll20ll . but few has been done to the massive Skyrme 
model, except for the exploration by Manton and Atiyah 111 711 . 

In this Letter, we revisit the static solutions of Symmetric 
Skyrmions in the S 1/(2) massive Skyrme model, which in- 
cludes the chiral symmetry breaking term (namely, pion mass 
term). An analtyical Skyrmion profile is proposed for the mas- 
sive Skyrme mode with the baryon number up to 8 by writing 
it as the hybrid form of a kink-like solution and the analytic 
solution obtained by the instanton method [6]. We show that 
this profile is in well consistent with the exact numeric solution 
obtained via the standard relaxation method. To determine the 
parameters involved in the profile, the downhill simplex algo- 
rithm is applied in lowering the energy as best as possible. The 
static properties of nucleon as well as delta are then computed 
using the B = 1 solution in the case of the reduced pion mass 



m = 0.526 within the framework of the semiclassical quan- 
tization of the Skyrme model, and compared to the previous 
computations and experimental data. 

2. The analytical profiles for Skyrmions 

The S U (2) Skyrme model [ 1 ] with pion mass term is given 



by 



L = J d 3 xi^tr(d ll Ud> 1 U i ) + -^tr [d M UU f ,d v UU f f 



iit„ i 
+-ffitr(U - 1) 

(1) 

with U (x, t) e S U(2) the nonlinear realization of the chiral 
field describing the cr field and n mesons under the constrain 
WU = 1, 2f„ is the pion decay constant, and e is the dimen- 
sionless constant characterizing nonlinear coupling. The last 
(mass) term in (Q} corresponds to the chiral symmetry breaking 
and ignored in the chirally symmetric Skyrme model. With the 
weak-7r expansion 

U = (cr + if. it)/f„ *1 + Jrr- 
Jn 

(t are the three Pauli matrices) the model ([TJ becomes the dy- 
namics of pian 



1 m 2 



in which the pion-mass term arises from the mass (third) term 
in (Q]). The Hamiltonian with respect to the Lagrangian ((TJ is 
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with m = nix/iefn) the reduced pian mass (in the unit of ef K ), 
and & are the partial derivatives in the length unit of 1 /ef„. 

We use the parameterization of the rational-map ansatz de- 
scribed in lfl5ll . 



U(r,z) = exp [iF (r)h R(i) ■ f] 



(3) 



where represents a map from the Riemann spherical coor- 
dinate z = e v tan((9/2), which is related to the angular coordi- 
nate (8, <p) via the stereographic projection, to the unit iso-triplet 
through the next stereographic projection in the internal space, 



1 



«R(z) 



1 + \R\ 



(2Refl, 2ImR, 1 - \R\ 2 ) . 



(4) 



Here, the map R(z) from a complex plane to S 2 = C U {00 j 
in the internal space has also to be a holomorphic map of de- 
gree B. This implies that R(z) can be written as a ratio of two 
polynomials R(z) = p(z)/q(z), where p and q are polynomials 
in z such that max{deg{p), deg(q)} = B, and p and q have no 
common factors. In 01, F(r) is a real profile function of the 
radial coordinate r, which has to satisfy the boundary condition 
F(Q) = n, F(oo) = 0. When 01 and <j4j applied, U becomes 



U(r,z) 



1 



1 + \R\ 2 



e iF + \R\ 2 e- iF 



2iR* sin(F) 
2iR sin(F) e~ iF + \Rfe 



2JF 



With 01 and (0, one has for the static energy for (0 



1 r 

37T Jo 



M 



dxx 
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Fl + 



;Bsm 2 (F)(l+F 2 ) 
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sin 4 (F) 



+ 2m 2 (l -cosF) 



(5) 



(6) 



(7) 



where 



4n , 



1 1 + kl 2 


dR 


^ 4 lidzdz* 


\l + \R\ 2 
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I d + kl 2 ) 2 



Here, the dimensionless variable x - ef K r and F x 
are used. The equation of motion of (0 is 



1 + 2B- 



sin z F 



2 sin(2F) 
Fxx + -F x + 



X \Fi-l- 1- 



_sin 2 F 



(8) 

dF(x)/dx 
(9) 



= m sin F, 



The numerical solutions to (|9]l and to the equation of mo- 
tion of the original functional (0 via the full field U calcula- 
tion are both presented by Battye, Sutcliffe and others [3]. It is 
shown there that an excellent agreement was achieved between 
the solutions given by rational map approximation and full field 
calculation for the baryon number B up to 12. Here, we re- 
examine only the rational map solutions to the Skyrmions with 
B = 1 ~ 8 for the case of the reduced pion mass m = 0.526, 
which is of the physically interest [ 10], by solving (0 analyti- 
cally and numerically. 



For given B, the minimization of the energy should be 
carried out in two steps: to minimize I (0 as a functional of 
R(z) at first and then to minimize energy for the optimized 
value of I. The first step is done in 112 111 and we simply list 
optimized values of I in Table I as a input for our calculation. 

The Skyrmion profile ansatz we propose is 



F(x) = 4w arctan [exp(-cx)] + 7r(l — w) 

,1/2 

^ / sinn i ///// v s 



X 



sinh (dmx) 



a 2 + sinh (dmx) 



(10) 



where m + 0, and (a, c, d, w) are the parameters to be deter- 
mined variationally. We note that the the second term in ( TTUb 
appears as the hyperbolic Skyrmion form obtained using instan- 
ton holonomvl lq 17 1 where a is related to the instanton scale 
A through 

'-<&>■ 

Given F(x) specified by (fTOt . we are able to minimize (0 
with respect to the variational parameters (a, c, d, w) for B = 
1 , ■ • ■ ,8 using the downhill simplex method (the Neilder-Mead 
algorithm). The optimized parameters are listed in Table I, in- 
cluding the soliton energy per B calculated accordingly and the 
corresponding numeric results obtained via the relaxation pro- 
cedure. 

The profile solution dTOb . with the parameters in Table I, is 
plotted in Fig.l and compared to the numerical solution to the 
equation (0. To the later, we employ the relaxation algorithm 
with the different initial setups of F(x), ranging from the Gaus- 
sian distribution, to the kink configuration, and find that the 
same result is achieved. It can be seen from Table I and Fig.l 
that a quite well agreement is achieved between the numerical 
solutions and that given by ( TTOb . 

We also present the fit of the parameters depending upon B, 
which are given by 

A = 0.4639 + 1.488 exp(15.038 - 15.3615) 
-0.5466 exp (1.7442- 1.78 16B), 



0.7196 + 2.287 lf^gf exp 



d = 0.7981fir u - 09 - 643; J 
w = 251.4- 



/6.774 



24X> 



l+exp{2 1.394 



fi+0.066519" 



20.404 



(12) 

and shown in Fig. 3. 

For the dependence of the charge one Skyrmion profile on 
the reduced pion mass m, we examine dTOb for the different 
m = 1,2,3, respectively. In Table II, the list is given for the 
parameters in HOi , which are calculated with the same proce- 
dure as before, and for the corresponding soliton energies ob- 
taind by (fTOlifmarked by Ana) and the energies obtained nu- 
merically (marked by Num). The chiral angle profiles ob- 
tained analytically and numerically are both plotted in Fig. 2 
for m = 0, 1,2,3, where the Skyrmion profile for the mass- 
less model is otained in [20] and used simply here for m — 
case. The results here show that the function dTOb for B > 2 
can be taken to be the approximated analytical solution to the 



2 



Skyrme model, which is presented in an analytical form for any 
low B. Since dTOb are not the most general expression of F(x), 
the variation of the parameters does not necessarily lead to a 
real analytical solution. However, considering that exact ana- 
lytical solution is difficult (almost impossible) to be found in 
most cases of the highly nonlinear problems like the Skyrme 
model, the solutions presented here can be an alternative to the 
exact analytical solution in so far as we are going to search the 
approximated solution which approaches the real solution as 
better as possible. 

3. The static properties of nucleons 

It is known that the static properties of nucleons can be ex- 
tracted by semi-classically quantizing the spinning modes of 
Skyrme Lagrangian using the collective variablesll9L llOtl . Fol- 
lowing Adkin et al. Ill Oil , we use the charge one solution of 
( ITOl ) to compute the static properties of nucleons and nucleon- 
isobar(A) within the framework of the bosonic quantization of 
Skyrme model. 

Choosing a S t/(2)-variable A(t) as the collective variables, 
and substituting U = A(t)Uo(x)A(ty into (Q]i, one can show that 
the Hamiltonian for the bosonically quantizatized model has an 
eigenvalue (H) = M SK + J(J + l)/(2/ A) with M SK the soliton 
energy for the static hedgehog Skyrmion(n«( Z ) = r) for which 
I = 1, 7 = /r/(3e 3 AX and 

Xoo 
x 2 dx sin 2 F [ 1 + F 2 X + sin 2 F/x 2 ] . (13) 

Hence, the masses of the nucleon and A-isobar are given by 
(with/ = 1/2 and 3/2) 

3 



M N = M* 



M A = M SK + 



8/oA 

15 
87oA" 



(14) 



As done in ifioll . we choose to adjust model parameters {f K , e) to 
fit the hadron masses, namely, the N, A , and n masses through 
(1141 1. where the other quantities involved are given by (0 and 
(1131 . Our reports for the computation of the static properties of 
baryons are presented in Table III, including the experimental 
data as well as the results in Ref.|9jjl0fl. The computations par- 
allel the computation reported in I1 1011 . We also list the most rel- 
evant formulas for our computations, for instance, the isoscalor 
root mean square(r.m.s) radius and isoscalor magnetic r.m.s ra- 
dius 



x 2 sin 2 FF X 



<W>B-{-§I 

1 J Q x 4 sin 2 FF x dx ^ 



M W ° \^x 2 sm 2 FF x dx\ 
and the magnetic moments for proton and neutron 



/=() , ../=! _ ( r >/=0 



Pp,n — P-p, n + Pp,n 



M N (M A - M N ) : 



M N 



2{M A -M N Y 
(15) 



where plus and minus correspond to proton and neutron, re- 
spectively. Since the chiral symmetry is broken explicitly, the 
isovector electric and magnetic charge radii are finite and can 
be calculated by 



e 2 /,V>/=i = | x 2 p,=i(x)dx 



Jr»oo 
x : 




(16) 



PMj=i(x)dx 



where pi=i(x) is the normalized electric iso-vector charge dis- 
tributions and can be derived from the zero component of the 
vector current. It can be shown that these two densities are 
identically same and can be given by 



pi=i{x) =Pm,/=iW 



x 2 sin 2 (F) {l + F 2 + sin 2 (F)/x 2 ) 
JT°° x 2 sin 2 (F) {l + F 2 X + sin 2 (F)/jc 2 } dx ' 



Therefore, the r.m.s radii of the isovector electric and magnetic 
charge are happen to be same. 

For the nuclear interaction, the pian-nucleon sigma term can 
be given by the matrix element in the nucleon states N{p) lfl8ll . 
which in our case becomes (see ifioll ) 



dV[2 - trU ] 



x dx( l - cos(F)) , 



where Uq is the the soliton solution of the chiral symmetric the- 
ory(namely, the massless theory). The result is cr — 38MeV, 
as given in 111 Oil . The calculated result for sigma term here is 
found to be 38.099MeV and listed in Table III, for which U is 
the the soliton solution of the chirally symmetry-breaking the- 
ory (namely, the massive theory ([TJ)- Being connected with the 
amplitude for meson-nucleon scattering, the sigma term can be 
extracted from the experimental data, with the "world average" 
estimate 36 + 20MeV 1181 [l9[], and is also listed in the table. 

4. Concluding remarks 

We show that the hybrid form of a kink-like solution and the 
approximated Skyrmion profile given by the instanton method 
are suited to approximate the exact profile for the symmetric 
Skyrmions in the Skyrme model with pion mass term for all 
baryon numbers up B to 8. This is first unified expression for the 
symmetric Skyrmions for the SU(2) Skyrme model with chiral 
symmetry broken. The computation of the corresponding soli- 
ton energies, and the static properties of nucleons within the 
framework of collective quantization, shows that the yielded 
Skyrmion profiles here are all in a good agreement with that 
given by the exact numeric solution. 
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The radial distance xi^ef^r) 
Fig. l.The chiral angle profile F (x) v.s. the radial coordinate for 
baiyon numbers from 1 (the most bottom) to 8 (the most top) and 
for the pion mass m=0.526. 



analytical 
numerical 



The radial distance x(=ef K r) 
Fig. 2. The chiral angle profile F(x) v.s. the radial coordinate for 
the pion mass from (the most top) to 3 (the most bottom) and 
for baryon number B= 1 . 



Table I 



B 


1 


2 


3 


4 


5 


6 


7 


8 


I 


1 


5.81 


13.58 


20.65 


35.75 


50.76 


60.87 


85.63 


A 


1.016 


0.3751 


0.4508 


0.4558 


0.4629 


0.4651 


0.4650 


0.4658 


c 


0.8058 


1.314 


1.306 


1.223 


1.166 


1.114 


1.070 


1.019 


d 


0.7981 


1.309 


1.237 


1.147 


1.093 


1.043 


1.001 


0.9532 


w 


3.203 


3.760 


20.45 


34.08 


58.63 


77.72 


85.91 


102.3 


|(Ana.) 


1.3084 


1.2848 


1.2603 


1.2078 


1.2230 


1.2147 


1.1817 


1.1984 


|(Num.) 


1.3077 


1.2833 


1.2593 


1.2069 


1.2219 


1.2132 


1.1799 


1.1953 



5 



Fig. 3. The values of the parameters A,c 3 d,w (the solid circles) 
presented in Tabic I plotted against B.Thc solid lines shows 
the fit of the corresponding parameters. 



Table II 



B=1,I=1 


m 


1 


2 


3 


A 


1.001 


0.9990 


0.9995 


c 


1.581 


1.813 


2.116 


d 


0.8814 


0.5419 


0.4236 


w 


-0.6877 


-1.746 


-2.209 


E(Ana.) 


1.418 


1.647 


1.855 


E(Num.) 


1.416 


1.642 


1.847 



Table III 



Quantities 


Ref.[9] 


Ref.flOl 


Analytical 


Num. 


Expt. 


M N (MeW) 


938.9(input) 


938.9(input) 


938.9(input) 


938.9(input) 


938.9 


M A (MeV) 


1232(input) 


1232(input) 


1232(input) 


1232(input) 


1232 


m^(MeV) 


O(input) 


138(input) 


138(input) 


138(input) 


138 


2A(MeV) 


129 


108 


108.20 


108.27 


186 


e 


5.45 


4.84 


4.8424 


4.8430 




<r 2 )ji 2 (fm) 


0.59 


0.68 


0.6825 


0.6821 


0.72 


<r 2 )^(fm) 


oo 


1.04 


1.0589 


1.0408 


0.88 


<'- 2 >M.t«,( fm ) 


0.92 


0.95 


0.9619 


0.9542 


0.81 




oo 


1.04 


1.0589 


1.0408 


0.80 


Hp 


1.87 


1.97 


1.9669 


1.9664 


2.79 




-1.31 


-1.24 


-1.2365 


-1.2369 


-1.91 


\Vp/Mn\ 


1.43 


1.59 


1.5906 


1.5898 


1.46 


o-(MeV) 




38 


38.099 


37.726 


36 + 20 
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